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Abstract 

In this paper we present a full set of 2- and 3-point functions for massless QCD at three-loop 
order in the MS scheme. The vertex functions are evaluated at the asymmetric point with 
one vanishing momentum. These results are used to relate the MS coupling constant to that 
of various momentum subtraction (mom) renormalization schemes at three-loop order. With 
the help of the known four-loop MS /3-function we then determine the four-loop coefficients of 
the corresponding MOM /3-functions. 

As an application we consider the momentum dependence (running) of the three-gluon 
asymmetrical vertex recently computed within the lattice approach by Ph. Boucaud et al. in 
[1]. An account of the four- loop term in the corresponding /3-function leads to a significant 
(around 30%) decrease of the value of the non-perturbative correction to the running 
found in [1]. 
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1 Introduction 



Momentum subtraction schemes provide a possibility to define a renormalization description for 

QCD in a rcgularization independent way. The concept of momentum subtraction is very old. It 
played an important role in the discussion of renormalization description dependence of physical 
quantities long ago [2, 3, 4, 5]. Recently, the momentum subtraction approach has been heavily 
used to relate lattice results for quark masses and coupling constants to their perturbatively 
determined MS counterparts (see, e.g. [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]). In [1, 13, 14] it 
has been argued that the knowledge of three-loop coefficients for the corresponding /3-functions is 
necessary and even the four- loop contributions should be taken in to account. This is because the 
accessible energy ranges in these calculations are just reaching a level where perturbative QCD 
calculations start to be valid approximations. 

A large subclass (in fact infinitely many) of mom schemes can be defined by subtracting vertices 
at the asymmetric point where one external momentum vanishes. We will call this point the zero 
point (ZP). These schemes are referred to as mom schemes and were introduced and discussed in 
some detail at two- loop order in [3]. A crucial fact for mom schemes is that setting one external 
momentum to zero for all three 3- vertices of massless QCD never will produce infrared divergencies. 

In this paper we present the perturbative calculation of the gluon, ghost and quark self-energies 
and all fundamental"'^ 3-vertices with one vanishing external momentum for massless QCD in a 
general covariant gauge. The one-loop triple gluon vertex was obtained in [17] in Feynman gauge, 
the result for general gauge can be found in [18]. At two loops, the triple gluon and ghost gluon 
vertices at the ZP have been determined in general gauge in [19]. References to earlier relevant 
publications and results for various momentum configurations can also be found in this work. The 
quark-gluon vertex can be found to two-loop order in Feynman gauge in [3]. 

These three-loop results at hand allow one to relate the coupling constants of any MOM-like 
scheme to the MS scheme at three-loop order. Recently, in [20] even the four-loop term of the MS 
/J-function has been computed. Using this result, we can also determine the /3-functions of any 
such MOM-schcmc up to (and including) four loops. 

The paper is organized as follows: In Section 2 the definitions of the ghost and quark self- 
energy and the gluon polarization are introduced and a brief outline of their calculation is given. 
In Section 3 we introduce the triple gluon, quark gluon and ghost ghion vertices at the ZP and our 
notation for these along with a description of their calculation. The calculation of the triple gluon 
vertex is performed directly and additionally in an independent way using the Ward-Slavnov- 
Taylor (WST) identity relating the triple gluon vertex to the ghost gluon vertex. This is content 
of Section 4. In Section 5 the MS renormalization procedure is described. In Sections 6 and 7 
we use the results of the previous Sections to obtain the coupling constants and /3-functions for 4 
different mom schemes of particular interest from their MS counterparts. Finally, in Section 8 we 
discuss the numerical importance of our results on the example of a recent lattice computation [1] 
of the momentum dependence (running) of the three-gluon asymmetrical vertex. 

The complete results are given in the appendix and also will be made available as input ffle 
for the algebraic programs form and mathematica at: 

http : // www-ttp . physik . uni-karlsruhe . de/Progdata/ 

^That is appearing in the QCD Lagrangian 
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2 The Gluon Polarization and the Ghost and Quark Self- 
Energies 

The QCD Lagrangian with Uf massless quark flavors in the covariant gauge is: 

n f 

C = -lG»,G'^A..^ + i^vi/[l>].,V'/-^(a^^P' + 9V(5^"-5rV^P, (1) 

The quark fields tpf transform as the fundamental representation and the gluon fields A" as the 
adjoint representation of the gauge group SU(3). and /"'"^ are the generators of the fundamental 
and adjoint representation of the corresponding Lie algebra. The 77" are the ghost fields and is 
the gauge parameter (^^ = corresponds to Landau gauge). 

From this Lagrangian one can derive three types of 2-point functions 



= DfM = -^j Axe^''^T[A1{x)Alm). 



where as usual r[Ai?] is the time ordered product of A and B and from now on we will skip 
the flavor index of the quarks. The propagators in eq. (2) can be expressed in terms of the 
corresponding self-energies in the following way: 



gab 

D""^ (a) = — 

t^u\q) (_^2) 



Sij (q) : 



(-q2)l + Sy(<z2) 



The self-energies Il{q^), ll{q^) and Sy(g2) can be calculated by applying the following projections 
to all one particle irreducible (IPI) diagrams with two external legs of the corresponding type: 





These projectors take into account that we are using dimensional regularization with D = 4 — 2e 
and are valid for a general SU(A^) gauge group. The trace has to be taken over Dirac matrices. The 
resulting scalar integrals are of the massless propagator type and can be evaluated with standard 
methods. For details about the technical setup used to perform these calculations we refer to the 
appendix. 
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3 The Triple Gluon, Ghost and Quark Vertex Functions 



In massless QCD there are the following 3-point functions: 

G(fj;'^(p,«) = Jdxdye-^^P-+'iy\T[A;{x)At{y)A;{0)]), 
G';!^"'%p,q) = j dx&ye-'^^^+'iy\T[ri'{x)n'^{y)Aim), 
Gfi-{p,q) = J dxdye-'^P-+''y\T[^i{x)^j{y)A;{0)]}. 

These are related to the vertex functions by 

G(1"'"=(P, q) = D^l, {-p)Dll, {-q)D;l, i-p q)r%^, {p, q, -p-q), 



Gf'''%P,ci) 

■,(3) a 



^"^-pK^P, q; -P - qW{q)DllAq + p), 



Gpij {P,q) = S,,'{-p)h\,,,,y{p,q-,-q-p)Sj'j{q)D'^^,^{p + q). 



(5) 



(6) 



Setting one external momentum to zero leaves only one momentum which can be used to construct 
tensor decompositions of the vertex functions. The color and Lorentz structure of all vertices at 
the ZP and our conventions are discussed in the next subsections. 



3.1 The Triple Gluon Vertex 

For symmetry reasons setting any of the three external gluon momenta of the triple gluon vertex 
to zero will lead to the same scalar functions. Up to two loops, the triple gluon vertex is known 
to be proportional to the totally antisymmetric color structure functions A color structure 

proportional to the totally symmetric cP^^ will not have a divergent part. In this work we only 
consider the /"''^ part. 

The triple gluon vertex needs to be totally symmetric under exchange of any two of the 
(bosonic) gluons and we are interested in the part proportional to f°-^'^_ So the Lorentz struc- 
ture for this part in the case of one vanishing external momentum is limited to the following 
three tensor structures, which are all antisymmetric under exchange of the two gluons with non 
vanishing momentum: 



r^t;(9,-9,o) 



'^gf {'^Qp.uqp - 9p.pq^ - 9puqp)Ti{q^) 



)qpT2{q^) + qp.quqpTsiq'^) 



(7) 



Due to the WST identity for the triple gluon vertex(for details see Section 4) the third function 
needs to vanish and finding T3(g2) = is a check for the correctness of the result. The functions 
Tj(g'2) can directly be calculated by applying the following projectors to the IPI diagrams with 3 
external gluon legs of which one carries zero momentum: 



Pit^^piQ) = 

P2llpi<l) = 

PSllpiQ) = 



abc 



7V(iV2 - 1) 

^ jabc 
• jabc 

N{N^ - 1) 



1 



D-1 
1 

D - 1 

qp,qvqp 



qiiq^qp Qupqp, gp.pqu 



q" 

iqp,q,^qp 

q' 



2^2 



2^2 



qp_ _ q^ _ qp_ 

9t^i' ^2 9t^P ^2 9i^P ^2 



ry6 



(8) 
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(9) 



The explicit calculation shows that indeed Ts vanishes. 



3.2 The Ghost Gluon Vertex 

The ghost gluon vertex is at tree level proportional to the momentum of the outgoing ghost. Since 
this vertex is the only interaction of the ghost field, it is clear that the vertex is also proportional 
to this momentum at any order in perturbation theory. This leaves two possibilities of one zero 
external momentum and due to the simple Lorentz structure also only two scalar functions to be 
determined (again we only consider only the /"'"^ color structure): 

ff^{-q,0-q) = -igr'^qM<l\ 

ff^{-q,q-0) = -igrt^fjq^). (IQ) 

The subscripts g and h stand for the external line that carries zero momentum; g for the gluon 
and h for the ghost. Again the Ti{q^) can be directly computed from the IPI diagrams with two 
external ghosts and one external gluon: 




rg(<?') 



be 



N{N^ - 1) g2 




(12) 



3.3 The Quark Gluon Vertex 

Finally the quark gluon vertex is proportional to the color structure T^" and for one vanishing 
external momentum there are two different possibilities (here setting either of the quark momenta 
to zero gives the same scalar functions up to three loops). A useful tensor decomposition is: 



A;„(-g,0;g) = gTtj 

A;,^{-q,q;o) - gi;;;. 



li^Agiq ) + 



A^(g^) 



(13) 
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where the additional subscript q corresponds to the case of a vanishing external quark momentum 
and T marks the transversal part. The A.^ can directly be calculated in the following way: 






-Tr 



NCf' 



1 



A{D-l) 



7m 
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(15) 



4 The Ward- Slavnov- Taylor Identity 

By using the WST identity for the triple gluon vertex, one can determine the scalar functions Ti 
and T2 from the gluon and ghost self-energies and functions related to the ghost gluon vertex. 
This check has also been performed in other determinations of the gluon vertex [17] [18] [19]. The 
general WST identity was determined in [21] [22] and can be found in the following form for the 
triple gluon vertex in [17] [18]: 



+J{q')G{e){gPq' - qPq,)fp^{q,k-p). 



(16) 



Here T^^j^{p, q; k) is related to the proper ghost gluon vertex (p, q; k) by: 

Tf^{p,q-k)^-igr''p''T,^{p,q-k). 
and we have introduced the functions 



(17) 



J(p2) = l + n(p2) and G(p^) 



l + n(p2) 
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In ref. [17] the following tensor decomposition for r^,^(p, q; k) is given: 

c{q,k,p) 

+%P,^d{q, k,p) +p^p^e{q, k,p). (18) 

The analytic determination of the full momentum dependence of the generalized ghost gluon vertex 
at three loops is impossible with the current calculation techniques. But of course we do not 
need to know the full momentum dependence for checking the WST identity at the ZP. Still, even 
the determination of all possible tensor structures for the case of one vanishing momentum and the 
corresponding expansions to first order in the vanishing momentum would be a very demanding 
project. To avoid the calculation of unnecessary parts of the ghost gluon vertex we follow closely 
the approach of [19] . 

For the vertex function there arc two independent momenta. They can be chosen in such a way 
that at the ZP one of them vanishes. Contracting the vertex in the WST identity (16) with the 
momentum that does not vanish at the ZP is relatively simple. In this case one can safely set the 
other momentum to zero on the right hand side and is left with functions of just one momentum. 
Inserting the tensor decompositions (7) and (18) gives: 

Ti(/)=a3(p2)G(p2)J(/) , as{p^) = a{0,p,-p). (19) 

There is also the possibility to contract the vertex with the momentum that vanishes in the limit 
of the ZP. This case gives a differential WST identity at the ZP. If the right hand side is expanded 
to first order in kp the constant lowest order term cancels. Dividing by kp, setting A: to on 
both sides and taking into account that for massless quarks G{0) = 1^ results in the following 
representation of this differential identity: 

^ij.upip, -P, 0) = +02(p^) J(/) {2gp,^pp - gp,pp^ - g^pPp) 
+^2^ )P [-^^^ + ^^"^^^ ) 

+d2(/)/j(/) (-2^^+g,pP. + g.pP,^ 

+J{p') (43(P') - 4i(P^)) - 9..Pp) ■ (20) 

where we have expanded a{p, k, q) to first order in small k and use the following shortcuts: 

d{p,0,-p) = d2{p^) , a{p,0,-p) = a{-p,0,p) = a2{p'^), 

a{-p-k,k,p) = a2{p^) + ^a'^s{p^) + 0{k^), (21) 

a{p,k,-p-k) = a2{p^) + '^a'^,{p^)+0{e) . 

Contracting cq. (20) with p'^ once more and projecting out the structure Ti(p^) gives another 
possibility to reproduce Ti[p^). It can be used to produce the simple relation 

a2(/) - p'd2{p') = G{p^)a3{p^). (22) 

Using this relation and projecting out the structure proportional to T2(p^) we can also obtain 
T2(p^) in a completely independent way: 

T2(/) = 2Ti(p^)-2a(/)A(/j(p2)) 

+J{p') {a',,{p') - a',,ip')) . (23) 

^G(O) = 1 since all contributing diagrams are massless tadpoles which vanish in dimensional regularization. 
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Finally, applying the projector P3 of eq. (8) to the right hand side of eq. (20) should give an 
alternative representation of T3{p^) and it is indeed found to be zero from this identity, as has 
been mentioned in the last Section. 

The scalar functions 02, as, and a'21 can be calculated from IPI diagrams with the following 
operations: 



be 



1 



(24) 



02(9^ 



pabc 




(25) 



a3((Z^) = l + P;r(<z)x 




(26) 



9.fcP,t^(g)rW X 




(27) 



\ 




(28) 



The additional gluon line is to visualize that at the outgoing ghost vertex we do not contract 
with the external momentum, which leaves an additional open index u as shown in eq. (17). The 



shortcut 7^ stands for a first order Taylor expansion in k and = <?a/3gf-gf-- Both are to 
be applied to the integrand before any of the integrations are performed. These manipulations 
are necessary because the mincer package can only handle scalar products of internal and one 
external momentum for more than one loop. Since there are only the two independent external 
momenta p and fc, the first order terms in k (of a scalar function) have to be proportional to k ■ p. 
This guarantees that the above manipulations will extract the scalars ajj. 

Using this strategy we found complete agreement for the bare expression with the ones from 
the direct computation of the triple gluon vertex. 
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5 Renormalization 



For a generic renormalization scheme^ R, the following relations between bare and renormalized 
quantities hold in dimensional regularization: 



(29) 



Here h = as/(47r) = g1/{lQn'^), ii is the 't Hooft unit of mass which is the renormalization point 
in the MS scheme. 

The self-energies and vertex functions are renormalized as follows 



Zf{h^, M, e) [1 + h^,e)] ,B^^..^^,_^B^^3,^ ■ 

Zf{h^, e) [1 + h^,e)] ,b=^..^,,,^b=^3^. 



zf(/i^M,e)f^>^/l^e) ^ 



(30) 



The renormalization constants for the self-energies can be read of eqs.(2,3). Likewise, from 
eqs. (5,6) and the QCD Lagrangian the renormalization factors for the trilinear operators can be 
related to the renormalization of the coupling constant and the field renormalizations 



yR d-B. yR 

zi z^- 



(31) 



QCD is known to be a rcnormalizablc theory. This means that Z^ can be found which make the 
renormalized Greens functions eqs. (2,5) finite when taking the limit e ^ and at the same time 
fulfill the WST identities, respectively eq. (31), to any order in perturbation theory. 

Using all these formulas, renormalization in the MS-schcmc is straightforward. It requires the 
Zx = Z^^ {h, e) to only contain poles in e and thus to be of the following form: 



Z.(/i,e) = l + ^izW(/i) , Zi'\h) = l + Y,h^Zi' 



(32) 



i>0 



j>i 



The /^-dependence of the fields and the coupling h is then described in the usual way by the 
renormalization group equations: 



m^^Mm) = 


-eh + l3{h) 


m = - 


OO 


u2 d (1) 

dh^^ ' 






73 (/i) = 


■DC 

i=0 


- h ^ Z^^^ 




= l3{h)'q{n) 


73 (/i) = 


i=0 




V^^A^(a*) = 




l2{h) = 


OO 

i=0 





(33) 



^In the following all quantities without explicit superscript refer to the MS scheme, the superscript B marks the 
bare quantities 
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where the last equahties can be derived from eqs. (29), the D dimensional /3- function \—eh + (3{h)\ 
and the fact that the bare quantities must be independent of ijl. As the functions in eq. (30) (3 
and the anomalous dimensions are here defined in D dimensions but are finite for e ^ 0. 

The MS renormalization constants and anomalous dimensions are well known up to three loops 
[23, 24] and for the /?-function [20] and the quark field anomalous dimension 72 [7] even the four- 
loop terms are known. Of course up to three loops the Zi can also be determined independently 
from the bare results for the self-energies and vertex functions. The corresponding anomalous 
dimensions are given in the appendix. 

Performing this standard renormalization procedure one arrives at the MS renormalized ex- 
pressions for all self-energies and vertex functions. The expressions for these in the limit e ^ 
at the point = —fi^ are given in the appendix for generic color factors and a generic gauge 
parameter. The momentum dependence can be obtained from the anomalous dimensions. The 
fact that indeed all poles in e cancel and eq. (31) holds is another check for the two- loop finite 
part and at least for the pole terms of any of the three- loop results. 



6 Four Particular mom like Scheme Definitions 



Momentum subtraction schemes are defined by setting some of the 2 and 3-point functions to 
their tree values for a fixed configuration of the states of the external particles (that is momenta 
and polarization state) in a certain gauge. This fixes the field renormalization constants at the 
renormalization point /z: 

1 + n"™' (-m' ) = 1 = Zf' (^2^ e) [1 + ( , e)] , 
1 + n'^'5M(_^2)l ^ ^ ^ Zf^{ix'',e) \l + II^(-/x2, e)l , 



l-hE^°"(- 



y,e) 



1 



(34) 



For the renormalization of the coupling constant there are infinitely many possibilities to define 
a momentum subtraction renormalization scheme, even when considering only the ZP. Not only 
is there an ambiguity in which vertex to subtract, but also there is the freedom to use a certain 
linear combination of the scalar functions appearing in the gluon and quark vertices, which can 
be related to fixed polarization states of the external particles [3] . If one considers the generalized 
ghost gluon vertex as defined in eq. (17) there are even more possibilities. 

Each of these choices defines a different renormalization scheme just as a different choice of 
the arbitrary renormalization scale. In general it is not possible to fix more than one vertex to it's 
tree value at the ZP with = — /x^, since they all are related by Ward identities. 

"The" MOM scheme originally was defined in [3] as a scheme in which all three triple-vertices 
are subtracted at the ZP without violating the Ward identities. This can be achieved by choosing 
a special set of the appearing scalar functions at one-loop order. Of course this universality does 
not hold anymore when considering the two- and three-loop order. This is why even in the original 
publication a generalization of this scheme for two and more loops is given. It subtracts the ghost 
vertex and corresponds to the original definition at one-loop order. Here we also consider the 
other two schemes that coincide with the MOM scheme at one-loop order. They are defined by 
subtracting just the scalar functions that appear at tree level of the triple gluon (we will refer to it 
as the MOMg scheme) and quark gluon vertex (quark momentum set to zero, MOMq), respectively. 

Another MOM like renormalization scheme definition (MOMgg) that uses the gluon vertex lately 
has been used to relate lattice results of the triple gluon vertex [12] and the gluon propagator [14] 
to perturbative calculations. The MOMq scheme that subtracts the quark gluon vertex has also 
been used in lattice calculations [16], but is not used beyond one-loop order in this work. 



6.1 Subtracting the Ghost Gluon Vertex 

A generalization of the original mom scheme to more than one loop is given by subtracting the ghost 
gluon vertex with the moment of the incoming ghost set to zero and a longitudinal polarization 
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of the external gluon (the case of a transversally polarized gluon will vanish at the ZP). The 
renormalization condition is: 



(35) 



Using cqs. (29-31,34), wc can connect the coupling constant in this scheme to quantities calculable 
in the MS scheme through the following chain of equations: 



Xm) 



1 
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= <?r(M)^r(M) 



(i + nB(p2))^i + nB(p2) 



(1 + n"^(p2))Vi + n"s(p2) 



(36) 



Squaring this and inserting the MS expressions for H, 11 and we arrive at the following relation 
between the coupling in this scheme and the MS coupling, expressed as an expansion in h = h"^ 
and* a = : 



h + h'^ 



+ 



169 10 9, 3^o] 
^ - yn, + -a + 4^!] + 

1913 100 2 117, , 351 



76063 4 , ^ ^ 



27 

42074947 , , ^ 3 , 836: 



1719 549 2 81 3 



16 



8362 , 



757 2 

1000 o 
n-^. 



769387 



2320 , 6931 , 



162 
4893 



16 2. . 



8505 4635 2 405 3 

- 2547 4^ 



3 

1485 



22287 2 21141 3 



28 2, 38 2. 199903 . 

:? - iH...3 _ 117 . ,4 _ 60675 
32 ^-^^^ 16 ^ 
290371 



64 



315 . 70245 

(37) 



6.2 Subtracting the Quark Gluon Vertex 

Using again the renormalization condition that at one-loop level is identical to the original MOM 
scheme we have: 



(38) 



This subtracts just the structure proportional to 7^^ which is present at the tree level and corre- 
sponds to a longitudinally polarized gluon. Similar operations as above then lead to 



and inserting the MS results gives: 
*see Section 7 for details 



Ar(-M') 



(1 + st^^(-M'))\/i + n^'"(-/x2) 



(39) 



11 



^MOMq 



h + +■ 



169 
1913 



10 



9 



K2 



27 



-Uf 



100 _ 



77035 4 ^ ^ ^ 
81^3 



42735097 333 , ^ 3 
^^728 2-"^^^^^ + 4"^' 



, 1647, 549,2 81 o' 

+ iscsa - 54C3 + -j^a + + 



3 , ,2 37579 , 2320 , 6759 



757 ,2 



386759 16 o , , 28 
n/ + yri/C3eL + 



16 



81 

1000 . 20067 , 2511 , . 

-4 + ^C3a + ^C3C^ 

255 



28 2. 
■n/C3 + 



729 

8325 3195 2 

-CbCl J^Cs^L 



16 

43311 2 20925 3 
^^2~^^ + ~6^^^ ■ 



16 



2547 4' 
64 



225 
~8 

315 , 



38 2^ lyyyu 



199903 o 



9 ' 

. 3 117 4 35385 
"C3Cl 32'''^^-^ 8 — 

.4 66765^ 281059, 
C5 + —^^L 



<3 



64 



(40) 



6.3 Subtracting the Triple Gluon Vertex I 

Subtracting just the Lorcntz structure that is present at the tree level the gluon vertex gives the 
following renormalization condition 



(41) 



This subtracts the triple gluon vertex with the zero momentum gluon and one of the others being 
polarized longitudinal (in direction of their momenta) and the last being polarized transversal and 
parallel to the plane defined by the gluon momenta (for details see [3]). A similar derivation as 
for the ghost gluon vertex leads then to the following relation of the coupling in the MOMg scheme 
to the MS coupling: 



(T'MS(_. 2\Y 



(i + n"s(_^2)) 



3 • 



(42) 



This gives the following relations between the coupling constants of the two schemes: 



^MOMg 



= h + 



169 10 



nr H £; 

12 9^2^ 

7571 100 2 81 , 207, 795, 549,2 81,, 

- + -W^'f + T^^^^ - — + -T^' + -16^- + 16^^ 

84464417 1375 , 3 , ,2 72161 , 2320 , 

+ ^456 ^"/^'^^^ + 4"^^^^^ + ^"^'^^ + ^"^'^^ 

39197 , 187 ,2 6098639 16 2 . . 28 2 . 

197149 2 1000 3 87729^^ 6261 ^2 1371 ^3 117 

-TOQRi O6805 4155 ss"^ ''1'^ 

64 ''^^^ ~ 64 

267071 

+ ■ 



38261 4 , ^ , 



2 , 38 2 ^ 
_ ■«/C3 + y 

1371 



279861 26805 4155 2 885 3 315 4 3285 



^7071 83511 2 41445 ^3 2547 4 
~6^^^ + ~6^^^ + ^28"^^ + "eT^^, 



(43) 
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6.4 Subtracting the Triple Gluon Vertex II 

The MOMgg scheme used in [12] and [14] is defined by the following renormahzation conditions: 

j.MOMgg^_^2^ - irf "^■s8(-^2^ = 1. (44) 

This subtracts the triple gluon vertex with the zero momentum gluon being polarized transversal 
and parallel to the plane defined by the gluon momenta. The other two gluon are also polarized 
transversal but perpendicular to the vertex plane. Another possibility to understand this linear 
combination is to realize that it corresponds to subtracting the transversal part of the 3-point 
function, which is the only one to survive at the ZP (see also [12]): 

G(f);^^(p, -p, 0) = {g,. - P^)p,G(3)(p2) (45) 
Inserting eqs. (6,3,7) shows that indeed: 



G(3)«bc _„ 0) - f« - Mfl)<o (Ti(p2) - ^T2(p2)) 

'^^iti^^^KP^ P'Uj - yg^^^^ ^ ''^''^(l+n(2)(p2))2(i + n(2)(0))' 

From this follows the relation of the coupling in the MOMgg scheme to the MS coupling: 



(46) 



^MOMgg ^ ^MS K^X y M J^2^2 IM )) _ ^4^^ 

(i + n"=(-/x2))' 



Inserting T\, T2 and 11 gives 
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(48) 



Following these examples and using the results in the appendix, it should be possible to con- 
struct the relations to relate the coupling constant in any MOM renormahzation scheme to the MS 
one. 



7 The MOM ;5-functions 

The first two coefficients of the /3-function in any massless renormahzation scheme are scheme 
independent. When considering renormahzation schemes in which the coupling constant depends 
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on the generic covariant gauge parameter, this statement is not true, since the additional gauge 
dependence spoils these universality just as a mass parameter [25]. This means that we need not 
only to take into account the gauge parameter when writing down renormalization group equations, 
but also that there is a difference between the gauge parameter defined in the MS scheme and in 
any of the MOM schemes we are considering: 



CMS ^3 (^MOM 



(49) 



In Landau gauge we do not need to consider the axiditional fj, dependence introduced by the gauge 
parameter nor the difference between and So in Landau gauge the /3- function in any 

MOM like scheme can be obtained in the following simple way once the MOM coupling is expressed 
as a series in the MS coupling 



Ah" 



dM2 



QljMOM 



(50) 



where one has to invert the series h^°^^^"^) and insert this series on the right hand side to express 
the MOM /J-function as a series in The four-loop QCD /3-function in the Ms scheme was 

obtained in [20] and reads: 
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(51) 



Combining this result with the relations between the coupling constants, we can finally obtain 
the first four coefficients of the /3-functions in the above introduced four mom schemes, where we 
use the definition: 



3 

E- 

i=0 



(52) 



and in Landau gauge /3o and f3i are independent of the renormalization scheme. The three- and 
four-loop contribution in Landau gauge for the four schemes under consideration are: 
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^MOMq 



29559 ■ 
+^ 594C3 



+ nf 



7769 
12 



64 ■ 



989 8, 

1- -Cs 

54 9^^ 



(57) 



Pa 



+ - 



+ ni 



2795027 174207 . 734415 



16 4 
737837 6709 



324 



27 



Cs 



32 

9280 



C5 



-Cs 



+ nf 



487751 67939 , 1016935 , 
+ — ^Cs + —rji—Cs 



12 

800 16 



144 



(58) 



^MOMgg 



186747 1683 
+~64 



+ nf 



35473 65 , 

1 Cs 

96 6 ^ 



+ n1 



829 8/ 



+ 



9 

(59) 



^MOMgg 



+ - 



20783939 



nf 



+ n} 



+ 4 



128 

2410799 

_ 64 
1464379 

648 
3164 



1300563^ 900075^ 
Cs — Cs 



32 
1323259 

^ 144 
12058 



32 



908995 , 

Cs + ^T^Cs 



144 



27 



Cs 



64 



27 



7540 
27 ^ 





320" 







(60) 



8 Discussion 



In numerical form the coupling constant relations and the corresponding /3-functions for these four 
schemes in Landau gauge read: 



^MOM ^ /i + [+14.0833 - 1.11111 n/] + /i^ [+475.475 - 72.4546 n/ + 1.23457 n^] 

+ [+18652.4 - 4109.72 n/ + 209.401 - 1.37174 n^] , (61) 

/juJokq ^ /i + /i2 [+14.0833 - 1.11111 n/] + /i^ [+470.054 - 72.4546 n/ + 1.23457 n)] 

+ /i^ [+18332.4 - 4089.25 n/ + 209.401 - 1.37174 n^] , (62) 

/jMMg ^ ^ _^ ^2 [+14.0833 - 1.11111 n/] + [+469.196 - 71.7046 n/ + 1.23457 n^] 

+ /i^ [+18332 - 4036.86 n/ + 206.568 n} - 1 .37174 n}] , (63) 

f^fMugg ^ ^ ^ [+23.3333 - 2.44444 n/] 

+ [+850.231 - 152.066 n/ + 4.64198 n^] 

+ [+37119.7 -9377.02 n/ + 602.56 n^- 7.7915 4] , (64) 

pilSk ^ (/i^i™^^ [-11 +0.666667 n/] + (/i'™'^ 

+ (h"°"y [-3040.48 +625.387 n/- 19.3833 n^] 

+ [h"°'''y [-100541 + 24423.3 n/ - 1625.4 + 27.4926 n^] , (65) 
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^ j^;ji55-Mg^ [-11 +0.666667n/]+ [-102 + 12.6667n/] 

+ (h'^^^y [-2971.42 +612.951 n/ - 18.8833 n^] 

+ (hf^'^sy [-99327.8 + 23904.4 n/ - 1591.07 + 27.0482 n)] , (66) 

^lETMq ^ (^/i-^^q^^ [-11 +0.666667n/]+ [-102 +12.6667n/] 

+ (h^^^y [-2980.85 + 621. 773 n/- 19.3833 n^] 

+ (h"°"^y [-98539.5 + 24249 Uf - 1622.19 + 27.4926 n^] , (67) 

^iSo-Mgg ^ (/i^^ss)^ [-11 +0.666667 n/]+ (/i^"sgj^ [_102 + 12.6667 n/] 

+ (y^'^ssy [-2412.16 +356.488n/ + 14.2834 n^- 0.888889 4] 

+ (^/iMOMggj^ [-84353.8 + 20077.1 n/ - 1433.44 + 108.637 - 3.95062 n)] . 

(68) 

The corresponding numbers for the MS /3- function are (as an expansion in h = h"^): 

[-11 +0.666667n/] + [-102 + 12.6667n/] 
+ /i^ [-1428.5 + 279.611 n/- 6.01852 n^.] 

+ [-29243 +6946.29 n/- 405.089 nj - 1.499314] . (69) 

It has been argued [2] that momentum subtraction schemes will lead to couplings that are only 
shghtly different. Comparing the MOMgg scheme with the other three schemes it is easy to see that 
this of course depends a lot on the Unear combination used in the renormalization condition. On 
the other hand it is clearly seen that the couplings as well as the /3-functions of the three schemes 
that are identical at one-loop order will also give very close-by values for the coupling constant at 
higher orders. Also the MOMgg scheme is still much closer to the other mom schemes than to the 
MS one. 

Inserting the values 1, . . . , 6 for n/ shows that the coefficients for all the values of n/ are smaller 
for the MS scheme than in any of the MOM schemes. The difference between the MS scheme and 
the MOM schemes is only small for n/ = 6 and gets quite large for the small values of n/. 

Recently the momentum dependence (running) of the three-gluon asymmetrical vertex corre- 
sponding to the combination Tgg = Ti{p^) — |T2 was computed within the lattice approach in [1]. 
The resulting (nonperturbative!) behaviour of 

a,(p2) = 47r/i^"sg(p2) 
in the flavorless QCD has been found to be best described by an ansatz^ 

«.(/)= «r*(p') + ^, (70) 

with 

= 237 ± 3MeV, c = 0.63 ± 0.03GeV^ (71) 
^To be consistent to [1] we use the Euclidean metrics below. 
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Table 1: A collection of fitted parameters obtained by imposing different values for the ratio 63/62 as 
defined in the text. "Whole" refers to the whole energy window (2 GeV < n <10 GeV) and "> 3" to a 
momentum range above 3 GeV. 



Here = exp{—7Q/%%)Kg;^^^^. The above results have been obtained by working at three-loop 
level and employing the MOMgg scheme. The authors of [1] have also investigated the dependence of 
the result on the (then unknown) four-loop contribution to /3"°"sg_ Their findings are summarized 
in Table 1, which we have copied from [1] with adding one more row obtained with the help of the 
linear extrapolation and corresponding to the true value of the parameter 

|i = ^OM = 2.78284... (f™n,=0). 

Tabic 1 clearly demonstrates that taking into account the four-loop term in the /3"°"sg_fun(2tion 
leads to a significant decrease (around 30%) of the value of the non-perturbative l/p^ correction. 



9 Conclusions 

In this paper we have analytically computed the full set of the three-loops propagators and funda- 
mental three-linear vertexes with one of three external momenta set to zero in the massless QCD. 
The results have been used to find the NNNLO conversion factors transforming the MS coupling 
constant to the ones defined corresponding to a set of regularization independent renormalization 
schemes based on momentum subtractions (MOM-schemes). Then we have used the conversion 
factors to evaluate the four-loop [3 functions in these schemes. 

The newly computed corrections to the coupling constant running in MOM-schemes prove to be 
numerically significant. They should be taken into account when confronting the running obtained 
with the help of lattice simulations with the pQCD predictions. 
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A Technical Remarks 



To evaluate the self-energies and vertex function at three loop order we made heavy use of computer 

programs. The program QGRAF [26] was used to generate the diagrams. For the self-energies and 
the vertex functions at the ZP we arrive at scalar integrals of the massless propagator type after 
applying some projectors. An algorithm which allows one to analytically evaluate divergent as 
well as finite parts of such integrals was elaborated in [27]. It has been implemented in an efficient 
way in the MINCER [28] package written for the symbolic manipulation program FORM [29]. The 
huge number of diagrams in the three-loop calculations requires a complete automation of the 
whole procedure which should include also the calculation of the color factors. This has been 
implemented and is described in more detail along with other similar installations in [30]. This 
setup also gives the user tools to perform series expansions in small parameters, a feature that 
was used to calculate the naive first order expansion in a small external momentum for the ghost 
vertex which is necessary to check the WST identity as described in Section 4. We had to evaluate 
the following diagrams: 
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All calculations were done on workstations using the Alpha 21164 processors running at 600 MHz. 
The approximate runtimes are only rough estimates and show that the calculation of the triple 
gluon vertex was by far the most demanding part. This is not only because of the number of 
diagrams but also because the complicated vertex and propagator structures generate a huge 
number of intermediate terms. 

Wherever possible, we have cross-checked the following expressions with the known two- and 
three-loop results, the gluon, ghost and quark field anomalous dimensions can also be found in [31] 
and [32]. In most cases we found complete agreement. We only find different values than [3] for 
the two- loop contributions to Ag and A^ (their Ts and r4). For the two-loop triple gluon vertex 
we find complete agreement with [19], which contains a list of miss-prints in earlier publications, 
among them the gauge dependence of the one-loop T2 in [3]. 

The latex code for all results in this paper has been created automatically. The code has then 
been retransformed to form input files and cross-checked against the original expressions. We 
hope that by doing these checks and avoiding any hand editing of the formulas we have circumvent 
the common problem of miss-prints in otherwise correct results. The price is of course that the 
layout of some of the formulas is not as fancy as it could be. Also note that these expressions will 
be made available in the World- Wide- Web. 

B Conventions 

B.l Color Factors 

The following results are given for generic color factors and are valid for any semi-simple com- 
pact Lie group, where T^^ are the generators of the fundamental representation and f"'^" are the 
structure constants of the Lie algebra. 
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B.2 Momentum Dependence 

We give the Ms renormalized result for all self-energies and vertex functions at the point — —(J?'. 
The correct dependence for any of the self-energies and vertex functions T can be restored from 
the renormalization group equations, which are with our conventions: 

^?-^T{h,q,^x,^L)= (^l^^-^+(i^+^3^L-^^r{h,q,ii,^L)=7rT{h,q,ii,^L) (B.l) 



For massless QCD one cam write 



r = ro + 5^r,W(iJ), (b.2) 

and as a direct consequence of (B.l) one gets (for n > 0) 

r„=(-/?|^-73aA+7r) Vr (B.3) 

The equation can bo directly used to reconstruct the g-dcpcndcncc of the scalar functions 1 + 11, 
1 + n and 1 + Sy. For the vertex functions one needs to multiply the functions T^, Ti and by 
g = AiT^/h. The anomalous dimensions of any F is: 

(rig nh~ nq \ 
7r = - ( + + J (B.4) 

where ng,nh and Uq arc the mimber of external gluon, ghost and quark fields, respectively, of the 
diagrams contributing to F. The QCD /3-function and the field anomalous dimensions are given 
in part D of this appendix. Note that in our definitions for the anomalous dimensions a relative 
sign and a factor 2 is different from many other publications. 



21 



C Propagators and Vertices in massless QCD 
C.l The Gluon Self-Energy 
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C.2 The Ghost Self-Energy 
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C.3 The Quark Self-Energy 
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C.4 The Triple Gluon Vertex 
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